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A formula for the Josephson current between two superconductors with anisotropic pairing sym- 
metries is derived based on the mean-field theory of superconductivity. Zero-energy states formed 
at the junction interfaces is one of basic phenomena in anisotropic superconductor junctions. In the 
obtained formula, effects of the zero-energy states on the Josephson current are taken into account 
through the Andreev reflection coefficients of a quasiparticle. In low temperature regimes, the for- 
mula can describe an anomaly in the Josephson current which is a direct consequence of the exsitence 
of zero-energy states. It is possible to apply the formula to junctions consist of superconductors 
with spin-singlet Cooper pairs and those with spin-triplet Cooper pairs. 
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I. INTRODUCTION 

The discoveries of the high- Tc superconductorsB have 
stimulated an intensive research in this field. A sym- 
metry of a Cooper pair is an important information 
to understand the mechanism of high- Tc superconduc- 
tivity. The Josephson effect in anisotropic supercon- 
ductors has attracted considerable interst in recent 
years because high- Tc supercojaductors may have the 
dx2_y2-wscve pairing symmetryaa. So far, transport 
properties in various junctions of the d-wave super- 

perconductors, a sign of the pair potential depends on a 
direction of a quasiparticlels, motion. As a consequence, 
zero-energy states (ZES's)Ell are formed at the normal 
metal/ superconductor (NS) interface when the poten- 
tial barrier at the interface is large enough. The ZES's 
have bee|B-|[Sfien in the conductance spectra of tunnel 
junctionsL^la. It is known that the ZES's cause a low- 
temperature anomaly of the JosephsO|ti_c.urrent in SIS 
junctions of the d-wave superconductomuj. 

The anisotropic superconductivity itself has been an 
important topic in condensed matter physics since uncon- 
ventional superconductivity was found in hea]rVj-fe.rrnion 
materials such as, CeCuzSia, UBeig and UPtaEEatJ. In 
a recent study, the anisotropic supercond-Hptivity was re- 
ported in a layered perovskite Sr2Ru04E3. Some of in- 
teresting effects of the anisotropy in the pairing sym- 
metrji-.|Ott— Josephson current are revealed in previous 
workEjtJ'El. However, in order to study the contribu- 
tion of the ZES's to the Josephson current, we have to 
pay careful attention to a boundasj condition of a wave- 
function at the junction interfaceE3. Thus an expression 
of the Josephson current that describes the effects of the 
ZES's is desirable to study an aspect of transport proper- 
ties in anisotropic supercondutor junctions. So far, such 
formular for the Josephson current-^s obtained in SIS 
junctions of dj;2_y2 superconductorgl3. However there is 
no general fomula which can be applied to junctions of 
spin-triplet superconductors. 



In this paper, we derive a formula for the Joseph- 
son current in junctions of anisotropic superconductors 
with spin-singlet and spin-triplet Cooper pairs. The re- 
sults are an extension of the Furusaki^Tsukada formula 
for s-wave superconductor junctionsE^. Effects of the 
ZES's on the Josephson current is naturally taken into 
account is, the obtained formula through the Andreev 
reflectionEa coefficients (ARC's) of a quasiparticle. The 
low-tcmpcrature anomaly in the Josephson current is de- 
scribed by the dependence of the ARC's on temperatures. 
Throughout this paper, we take the units of ?i = fc^ = 1, 
where is the Boltzmann constant. 

This paper is organized as follows. In Sec. II, we derive 
the Josephson current formula based on the mean-field 
theory of superconductivity. In Sec. Ill, the formula is 
applied to junctions of superconductors with spin-singlet 
and spin-triplet Copper pairs. The conclusion is given in 
Sec. IV. 



II. JOSEPHSON CURRENT FORMULA I 

Let us consider SNS junctions as shown in Fig. ||, where 
the length of the normal metal is Tat and the cross section 
of the junction is Sj. The (BCS) Hamiltonian in the 
mean-field approximation reads 
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where Ccr{r) is the annihilation operator of an electron 
at r with spin a or |, {c(r)}' is the transpose of 



Eq. (pT^), (7o is the unit matrix of 2 x 2, and /xj? is the 
Fermi energy. Spin-independent potential is represented 
by Vo{r) which includes the barrier potential at the two 
NS interfaces given by Vf, {S{z) + S{z — Ln)}. Spin-orbit 
scattering in the normal metal is denoted by V{r) ■ &. A 
pair potential between an electron with {a, r) and that 
with (a', r') is described by Ao-^o-'li — i"'). In the normal 
segment (0 < z < L^), the pair potential is taken to be 
zero. In what follows, 2x2 matrices are indicated by — . 
The pair potential is given by 



A(r) 



ida{r)a2 singlet 
i{d{r) ■ cr)d-2 triplet 



(2.4) 



where dj with j = 1,2 and 3 are the Pauli's matrices. 
The pair potential satisfies a relation 



-A*(r'-r) = A(r-r') 



(2.5) 



The Hamiltonian in Eq.(2.1) is diagonalized by the Bo- 
goliubov transformation, 
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denotes the annihilation operator of a Bogoliubov quasi- 
particle. The wavefunctions satisfy the Bogoliubov-de 



Gennes (BdG) equationEl, 
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When the wavcfunction 
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is belonging to a positive eigenvalue Ex, the wavcfunction 

(2.12) 
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is belonging to —Ex- They satisfy the following relations 
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FIG. 1: The SNS junction of anisotropic superconductors is 
illustrated. The phase of the pair potential on the left (right) 
superconductor is (<^fl)- 



where X)a' ^ summation over A with positive eigenval- 
ues. The local charge density is defined by 



(2.17) 



where i is a time. The current conservation low implies. 



d 



dt 



-^P{r,t)+\I ■ J{r,t) ^ 0. 



(2.18) 



The Josephson current between the two superconductors 
is calculated from the expectation value of Eq. ( 2.18| ) 
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where T is a temperature, Q^^{r,r')jsthe Matsubara Green function of the SNS junctions and indicates 4x4 
matrices. On the derivation of Eq. (2.19), we have assumed that the amphtude of the pair potential is much smaUer 
than the Fermi energy ^p. 



In the superconductors, we assume that the all potentials are uniform. Thus the BdG equation in Eq. (2.10) is 
given in the Fourier representation. 
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where = /{2m) — iip and 
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Since relations 



are satisfied in the momentum space, one finds 
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When z < 2' < 0, the Green function can be calculated as 



X <i> 







{ 








( 









kp{k\,z)ai 



+ ]K,{-kt.z)+[ X,(-fcl,2)a2 ;t \K,{kl,z)b2\K,{kt,z') 



kl_ 



n: 



(2.27) 



with 



h 



Ah 



2m{^iF - e{p) + i^i,±), 
2m{nF - e(p) - i^i,±), 



f^;,±=xMU|Ai,±|2, 



'2m' 
do,± =do(p, ±^2), 

d± ^d{p,±kz), 



ic?Q -l-(T2 singlet 
i{d± ■ &)d'2 triplet 



(2.28) 
(2.29) 
(2.30) 

(2.31) 
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(2.33) 
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where Lpj for j = L or R is the phase of the superconduc- 
tor, p = {kx,ky) and p — {x,y). The amphtude of the 
pair potential for unitary states is defined by 



lA 



i.±\ 



IA4 
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\d±\ triplet 



(2.40) 



In unitary states, these amplitudes are independent of I, 
where I indicates the spin configuration of a quasiparticle. 
The amplitude of the pair potential depends on the spin 
configuration of a quasiparticle in nonunitary states, 



I A, 
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In Eqs. ( 2.28| ) and ( ^.29 ), ki^^'' is the wavenumber in 



the electron (hole) branch for l-th spin state. In the fol- 
lowing, we approximately describe these wavenumbers as 
fcf'^'' K kz = ^j2m{iiF — e(p)) as shown in Eqs. ( 2. 32 ) 



and (2.33), where (p, ifc^) is the wavenumber on the 
Fermi surface. The /-th column of 
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Throughout this paper, we use a representation 
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for unitary states. In unitary states, Ai^± is iadependent 
of I because of q = 0. For nonunitary stateslEl, we use 



u± 



(2.50) 
(2.51) 



corresponds to the wavefunction of ^-th spin state in the 
electron (hole) branch. The reflection coefficients from 
the left superconductor |-tf) the left superconductor are 
defined in a matrix formE3 



aj=i,2 - 
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The ARC from the Z-th spin state in the electron (hole) 
branch to the Z'-th spin state in the hole (electron) branch 

is denoted by ai{l' , I), {0,2(1', I))- In the same way, &i(/', I) 
(62 (^',0) is the normal reflection coefficient from the l- 
th spin state in the electron (hole) branch to the /'-th 
spin state in the electron (hole) branch. These reflection 
coefficients depend on p which indicates the propagating 
channel a t the left NS interface. Substituting Eq. (2.27) 
into Eq. (2.19), the Josephson current becomes 



ii'l 



(2.45) 



The expression of the Josephson current in Eq. ( 
is an extension of the Furusaki-Tsukada formuL 
s-wave junctions. 
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(2.52) 
(2.53) 
(2.54) 

(2.55) 

(2.56) 
(2.57) 
(2.58) 
(2.59) 
(2.60) 
(2.61) 



In this paper, we consider four reflection processes to 
calculate ai and 0,2 as shown in Fig. ^ (a) and neglect all 
higher-order terms. This approximation is justified when 
the potential barrier at the NS interfaces is large enough 
and the transmission probability in the normal segment 
is low enough. Thus in the normal segment, insulators or 
dirty normal metals are assumed. In order to estimate di 
and 02, we calculate the transmission and the reflection 
coefficients at the single NS interface for fixed p as shown 
in Appendix A. The ARC's in Fig. (a) are given by 
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a(i) (p) = is%iP, L) ■ il^, ■ f^V(p', R) ■ ip',p ■ tNsiP, L), (2.62) 
p' 

a'^^\p) = J2 is%iP, L) ■ il^p, ■ f%\{p', R) ■ . L), (2.63) 

p' 

ai'\p) - E ^""a^(P' ^) • ^P.P • ^^'a'(p'' ^) • ^P'.P • ^"^'^^(P' (2.64) 
p' 

(P) = E ^"I'a^(P' ^) • ^"p,P' • ^nn{p', R) ■ %',p ■ ifsiP, L), (2.65) 
p' 

where ip/^'p is the transmission coefficient of the electronlike (holelike) quasiparticle in the normal conductor, and p' 
indicates the propagating channel at the right NS interface. The transmission coefficients in the normal metal are 
described by 

t>,p =ivp e-^'=^^« I dp I dp'g^fip',LN;p,0)x;'ip')Xpip), (2-66) 
f^,p, =zty e^'^^^" J dp I dp'g^;''ip,0;p',LN)x;ip)Xp'ip'), (2-67) 



where gu^'^^^\r , r' ) is the Green functiorL-Qf the nor- 
mal conductor in the electron (hole) branchcZl. The ve- 
locity of a quasiparicle in the z direction is Vp for the 
propagating channel with p. We assume that the NS 
interface is sufficiently clean so that p and p' are con- 
served while the transmission and the reflection at the 
interfaces. In aj^'' in Eq. ( p.62| ), a quasiparticle- wave is 
initially incident into the normal segment from the left 
superconductor through the channel specified by p. Af- 
ter the Andreev reflection at the right NS interface, we 
assume that the reflected wave transmits to the left su- 
perconductor through the initial channel of p. This is 
because a quasiparicle in the normal segment has the 
retro property under the time reversal symmetrica. The 
two arc's in Eq.( 2.45 ) are given by ai = a!]^^ + a^^-* and 

Cl2 — CI2 ' ^2 • 



By using Eqs. (2.45) and (2.62)-(2.65), we can derive 



a general expression of the Josephson current 

P p' 
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The reflection processes in Eq. ( 2.68 ) are summarized in 
Fig. (b). Since the relations 
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are satisfled (see Appendices A and B), the Josephson 
current results in 



J==-2eImEE^E 
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The formula in Eq. (2.73) can be applied to various 
Josephson junctions. For instance, it is possible to cal- 
culate the Josephson current in clean SIS junctions by 
using a relation i'^'p, oc dp^piaa- We also note that the 
two superconductors are not necessary to be identical to 
each other. 



III. JOSEPHSON CURRENT FORMULA II 

In this section, we show the ARC's of the superconduc- 
tors in spin-singlet, spin-triplet unitary and spin-triplet 
nonunitary states because the Josephson current is de- 
scribed by the ARC's at the NS interfaces in Eq. ( |2.73| ). 

Firstly, we consider the superconductor with the spin- 
singlet Copper pairs. The ARC's are given by 
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FIG. 2: Four reflection processes in (a) contribute to the 
Josephson current. The Josephson current calculated from 
the four reflection processes in (a) is summarized in the re- 
flection processes in (b). 



where H = mVb/kp represents the strength of the poten- 
tial barrier at the NS interface and j = L or R symboh- 
cally denote the character of the superconductors such as 
symmetries of the pair potential and orientation angles. 

Secondly, the ARC's in spin-triplet unitary states are 
given by 



~eh 



' NN 



{p,L) = - ittu{p,L) e 
{p,R)=-irl{p,R) e-'^«, 



ftn(p, j) =iTtu{p,j) ■ a- &2, 
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(3.9) 
(3.10) 

(3.11) 

(3.12) 
(3.13) 



In unitary states, d± often has a single component. In 
such case, one finds 



Ttu{p,j)^klK^ 



d 



(3.14) 



because of dl x d_ = 0. 



Finally we show the ARC's in nonunitary states. 
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(3.21) 
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Detail of the calculation is shown in Appendix A, where 
we derive the ARC's of the superconductors in nonuni- 
tary states. We do not show the derivation for unitary 
states because it is much simpler than that in nonuni- 
tary states. As shown in above equations, the expression 
of the arc's in nonunitary states is very complicated. 
However if a relation 



d =d^ = i^d , 
ly =1 or — 1, 



(3.23) 
(3.24) 



is satisfied, the ARC's can be reduced to a rather simple 
expression 



s„„(z) ^Hmi-iy)\LOr,\ + {i + i^)ni} 

+ Pj\uJn\+ni). 



(3.25) 
(3.26) 

(3.27) 



The effects of th e Z ES 's on the A RC's can be easily con- 
firm ed in Eqs. (B.5), (3.12) and ( [3.27 ). For instance in 
Eq. (3.27), we find in the limit of i7 >> 1 and ujn 0, 



^nu (0 




(3.28) 



In the absence of the ZES's {i^ ~ 1), the reflection co- 
efficients proportional to On the other hand in 
the presence of the ZES's {v = —1), the reflection coeffi- 
cients are independent of the barrier height. In this way, 
the low-temperature anomaly of the Josephson current is 
described by the ARC's. 
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In the normal metal, the two Green functions in 



Eqs. ( 2.66 ) and ( 2.67 ) satisfy a relation as shown in Ap- 
pendix B, 



0-2, 



(3.29) 



because of the time reversal symmetry. The transmission 
coefficients can be parameterized by 

t{p',p) =To{p',p)ao +t{p',p) ■ a, (3.30) 

=^Avv / / '^p'x*p'{p')xpip) 

xg^;'{p',LN;p,0). (3.31) 

Since the amplitude of the spin-orbit scattering is much 
smaller than that of the spin-independent transmission 
probability, we assume that 



M » |T|. 



(3.32) 



The conductance of the normal metal at T = is given 
by 



Gn = lim J- Tr Vf(p',p)ft(p',p), (3.33) 
p,p' 

=^ lini^^5]|ro(p',p)p. (3.34) 



By using Eq.(3.30), the Josephson current is rewritten 



as 



J = -2eIm^^T^ Tr 

X [r'NN{P.L)-&2{Ta(ro + 



'NN 



{p',R) ■ {roo-Q 



(3.35) 



Firstly we consider Josephson junctions where the two 
superconductors have the spin-singlet Cooper pairs. The 
Josephson current is given by 



Jss =4e sin '^T^^ 

LOn p,p' 

xr,„(p',i?)|ro(p',p)pr,„(p,L), 



(3.36) 



where f = (fiL - 'PR- 

Secondly we consider junctions where spin-triplet and 
spin-singlet superconductors are on the left and on the 
right hand sides, respectively. The Josephson current 
results in 

Jts =4eT^^Im 

UJn P,p' 

X [ e^^r,„(p', R)W{p',p) ■ Tt{p, L)] , (3.37) 



W{p',P) = {toT + tqt* + IT* X t) (p',p) 



(3.38) 



where Tt represen ts Ttu in Eq. (3.11)orr„„ in Eq. (3.18). 
As shown in Eq. ( 3.38 ), the Jts vanishes in th e abs pace 
of the spin-orbit scattering in the normal metalEaESlti. 

Finally when the two superconductors have spin-triplet 
Cooper pairs, the Josephson current is given by 



Im 

t^Tl PiP' 

x[e^^Tt{p,L)-T:{p',R)\T^{p',p)\' 



(3.39) 



The obtained formula in Eqs. ( |33^ ), and 



are essentially the same as those in the previous resultsEil 
when the ZES's are not formed at the NS interfaces. 
However in the presence of the ZES's, the dependence 
of the Josephson current on temperatures in our results 
is drastically different from that in the previous one's. 
This is because the ARC's ( F^u, Tm and r„„) describe 
the low-temperature anomaly of the Josephson current 
in the SNS junctions of anisotropic superconductors. 



IV. CONCLUSION 

On the basis of the mean-field theory of the supercon- 
ductivity, we derive a formula for the Josephson current 
between two anisotropic superconductors. The Joseph- 
son current is expressed by the Andreev reflection coef- 
ficients at the junction interfaces. The contribution of 
the zero-energy bound states formed at the NS interfaces 
to the Josephson current is taken into account through 
these Andreev reflection coefficients. The formula can be 
applied to SIS and SNS junctions of the anisotropic su- 
perconductors with spin-singlet and spin-triplet Copper 
pairs. 
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APPENDIX A: TRANSMISSION AND 
REFLECTION COEFFICIENTS AT THE NS 
INTERFACE 

We derive the transmission and the reflection coeffi- 
cients at the left NS interface {z = 0), where the super- 
conductor is in spin-triplet nonunitary states as shown in 
Fig. ^. In what follows, we calculate the coefficients after 
the analytic continuation (i.e., E — > iLOn) for w„ > 0. In 
the normal metal, a wavefunction of a quasiparticle can 
be described by 
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-ikz z 



-ikz z 



Xp(p), 



(Al) 



where a and (3 {A and B ) are the ampHtudes of incoming (outgoing) waves in the electron and the hole branches, 
respectively. In the same way, a wavefunction in the superconductor is given by 



u+ 



Xp{p), 



(A2) 



where 7 and S (C and D ) are the amplitudes of incoming (outgoing) waves in the electron and the hole branches, 
respectively. We note that a, /3, 7 and 6 have only diagonal elements. 

The two wavefunctions satisfy a continuity-condition at the left NS interface. 
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(A4) 



From Eqs. (A3) and (A4), we obtain the transmission Here we define 
and the reflection coefhcients / 
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I u 
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(A6) 
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(A8) 
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(All) 
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Zl {u'_y^Cl^ e*^^/2, (A14) 



(A15) 

(A16) 
(A17) 
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l,± 



K —kz + iH. 



(A18) 

(A19) 

(A20) 

(A21) 
(A22) 



In the same way, the ARC's at the right NS interface are 
given by 



{-p,R)^[f%^AP,R)\ 
On the derivation, we use identities. 
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2Q± 

Pl.± 



Sli' 



2\Q±\Ql 



Sl^ • A± • Al • Si>,± 



Ai ■ Si,± ■ Sl ± • A± 



iiS. 
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^^^^ 
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(A23) 
(A24) 

(A25) 

(A26) 
(A27) 
(A28) 

(A29) 

(A30) 
(A31) 
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The arc's of superconductors in unitary states can be 
calculated in the same way. The derivation of the ARC's 
in unitary states is much simpler than that in nonunitary 
states. 



Superconductor 
9t 



Normal Conductor 





z = 

FIG. 3: Amplitudes of incoming and outgoing waves at the 
left NS interface. 

In addition to the four reflection processes shown in 
Fig. ^(a), six reflection processes can be considered for 
Si and 0,2 a s shown in Fig. ^ By using the coefficients 
in Eqs. ( |A5[ )-( A16), it is possible to show that these six 
processes do not contribute to the Josephson current. 
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APPENDIX B: TRANSMISSION COEFFICIENTS 
IN NORMAL METAL 



Since the amplitude of the pair potential in the normal 
metal is taken to be zero, the BdG equation in Eq.(2.1C) 
is decoupled into two equations, 



-hQ{r)vx =v\E\. 



(Bl) 
(B2) 



The Green function in the normal metal obeys the equa- 
tion, 

{iujnao - ho(r))g^f{r, r') ==,5(r - r')<To, (B3) 
{iu;„ao + h*{r))g^;'^{r, r') =5{r - r')a^. (B4) 

The Green function in the two branch are represented by 



4(r'), (B5) 



(B6) 



where we use the complex conjugate of Eq. (Bl) f or the 
Green function in the hole branch. By using Eqs. (2.66) 
and (2.67), we can show relations 



Wp'.I- 



(B7) 
(B8) 



When the time-reversal symmetry holds in the normal 
metal, we find 
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3(5) «'- - 
h 




' NN 


^NS 
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'NN 






^ eh 
hN 



FIG. 4: Reflection processes included in the coefficients ai 
and 5,2 ■ These processes, however, do not contribute to the 
Josephson current. 



hl{r)ia2U\ = ia2UxEx. 



(B9) 



The Green function in the hole branch is described by 
that in the electron branch. 



G5;ir',r) = -a2 gZf{r,r') a2 



(BIO) 
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